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Abstract 



We analyze the one loop corrections to hadronic Z decays in an R-parity 
violating extension to the Minimal Supersymmetric Standard Model (MSSM). 
Performing a global fit to all the hadronic observables at the Z-peak, we 
obtain stringent constraints on the R- violating coupling constants A' and A". 
The presence of these couplings worsens the agreement with the data relative 
to the Standard Model. The strongest constraints come from the b asymmetry 
parameters A^ and Ap^ib)- From a classical statistical analysis we find that 
the couplings A '3^, A '32, and A321 are ruled out at the la level, and that A^33 
and A33j are ruled out at the 2a level. We also obtain Bayesian confidence 
limits for the R- violating couplings. 
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I. INTRODUCTION 



R-parity conservation is often assumed in supersymmetric model building in order to 
prevent a host of phenomenological complications such as fast proton decay. This also 
serves to make the lightest supersymmetric particle (LSP) stable and thus provide a dark 
matter candidate. However, R-parity conservation is not a necessary condition for avoiding 
many of these problems. For example, the imposition of other discrete symmetries, such 
as conservation of either baryon number or lepton number, may be adequate to provide 
phenomenologically acceptable models. (For recent reviews, see Ref. [Q.) Furthermore, the 
evidence for neutrino mass recently observed at Super-Kamiokande lends improved mo- 
tivation to consider R-parity violating extensions to the minimal supersymmetric standard 
model (MSSM). Therefore, one is led to question just how much R-parity violation can be 
introduced without conflict with current experimental data. In this paper, we study the 
radiative corrections from R-parity violating extensions of the MSSM to the electroweak 
observables in hadronic Z decays, namely the ratios of hadronic partial widths and the par- 
ity violating asymmetries. Experimental data from LEP and SLD place stringent limits on 
the size of these corrections, thereby constraining the possible strengths of the R-violating 
interactions. 

We focus on the effects of the R-parity violating superpotential and neglect possible 
effects from the corresponding soft-breaking terms 0. This simplification allows us to 
rotate away the bilinear terms In this case, the R-parity violating superpotential has 
the following form: 

= 2^ijkLiLjEk + X'ijkLiQjDk + -\'-ji.Uif)jDk , (1.1) 

where Li, Ei, Qi, Ui, and Di are the MSSM superfields defined in the usual fashion [Q, 
and the subscripts i,j,k = 1,2,3 are the generation indices. These interactions can give 
potentially sizeable radiative corrections to the hadronic observables depending on the size 
of the coupling constants A, A', and A". The A couplings are already tightly constrained 
to be O(10~^) or less, and their effect on the Z-peak observables is negligible f^.Q The 
constraints on the A' and A" couplings are much less stringent. However, they cannot be 
present simultaneously in the Lagrangian since this would lead to unacceptably fast proton 
decay Therefore, we can make the further simplifying assumption that only one or other 
of the operators LiQjDk and UiDjDk is present at a time. 

When constraining R-violating interactions using experimental data, it is important to 
provide a consistent accounting of the corrections from the R-conserving sector also since 
they may be sizable depending on the choice of SUSY parameters. It is also important to 
include all the affected observables in a global fit since different observables may pull the 
fit values in opposite directions. This was illustrated in our previous paper in which 
the violation of lepton universality was used to constrain the A' couplings. There, the Z- 
lineshape observables alone preferred a 2a limit of jAgg^l < 0.30, but a global fit resulted 



They do not affect the quark couplings in any case. 
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in |A33^| < 0.42. Neither of these points were considered in previous works such as Ref. Q 
where R-conserving corrections were neglected altogether, and only corrections to the ratios 
of hadronic to leptonic partial widths = Fhad/ra = e, fi, r) were considered. It is clear 
that these ratios receive R-conserving corrections from top-Higgs and chargino-sfermion 
loops, as well as QCD and gluino corrections which depend strongly on as(M^). Therefore, 
the resulting la bound of jAg^^l < 0.50 of Ref. P] is hardly robust. 

In this paper, we consider all the purely hadronic observables which can be expressed as 
ratios of the quark couplings to the Z, i. e. the ratios of hadronic partial widths and the parity 
violating asymmetry parameters. These are unaffected by QCD and gluino corrections since 
they modify the left- and right-handed quark couplings multiplicatively, leaving the ratios of 
the couplings intact .0 The rest of the R-conserving sector induces relevant corrections to the 
left-handed quark couplings only, whereas the R-breaking sector affects predominantly the 
right-handed quaik couplings. This allows us to parametrize and constrain the R-conserving 
and R-breaking corrections separately, thereby constraining the R-breaking sector without 
making ad hoc assumptions about the R-conserving sector.^ Also, since we incorporate into 
our fit the corrections to the forward-backward and polarization asymmetries which are 
much more sensitive than to the shifts in the right-handed quark couplings, we are able 
to substantially improve the limits on the R-breaking interactions. 

This paper is organized as follows: In section |I| we discuss the approximations we make 
to simplify our analysis. In sections |T| and we discuss how the A' and A" interactions 
affect the couplings of the quarks to the Z. Section discusses the corrections from the R- 
conserving sector. In sections ^ and [VI I| , we parametrize the R-conserving and R-violating 
corrections to the LEP/SLD observables and fit them to the latest expermental data, and 
then translate the result into limits on A' and A". In section VIII , we provide the Bayesian 
confidence limits on A' and A" with the a priori assumption that the MSSM with R-violation 
is the correct underlying theory. Section 11^ concludes. 



II. PRELIMINARY SIMPLIFICATIONS 

As we stated in the introduction, we only consider supersymmetric R-violating inter- 
actions and neglect the effects of soft-breaking R-violating terms.0 We also neglect the A 
interactions and consider only the A' or the A" interactions at a time. In addition, left-right 
squark mixing is neglected since their effects are expected to be unimportant 0. Even with 
these simplifications, we still have 27 independent A' couplings or 9 independent A" couplings 
which must be considered. 



^We assume degenerate squark masses. 

^Similar methods have been used in Ref. and to constrain flavor specific vertex corrections 
while taking into account the flavor universal oblique corrections. 

^See Ref. H for a discussion on their possible effects. LEP/SLD hadronic observables can also be 



affected by resonant sneutrino production |12|. 
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However, a careful look at the diagrams which must be calculated begets a further 
simplification. The corrections to the Zqq vertex generated by the A' and A" interactions in 
the superpotential fall into four classes: 

1. One particle irreducible (IPI) diagrams with two scalars and one fermion in the loop. 

2. IPI diagrams with two fermions and one scalar in the loop. 

3. IPI diagrams with two fermions and one scalar in the loop, with two mass insertions 
on the fermion lines. 

4. Fermion wavefunction renormalization diagrams. 

In these diagrams, it is clear that the scalar must be the sparticle while the internal fermion 
must be an ordinary lepton or quark. The invariant masses of the external gauge boson and 
the external fermions must be set to m| and ~ 0, respectively. 

Of the four classes, the third class is finite while the 1/e poles of the first two classes cancel 
against the poles in the fermion wavefunction renormalizations. An explicit evaluation of the 
finite pieces of the diagrams reveal [0 that they lead to numerically significant contributions 
only when the fermion running in the loop is heavy. In fact, the amplitude of a diagram with 
a massless internal fermion is only about 10% of that with an internal top-quark, assuming 
that the scalar (sfermion) mass is the same. Since each diagram is proportional to A' or A" 
squared, dropping these 10% contributions to the amplitude will result in a 5% uncertainty 
in the limits obtained for the A' and A". We can therefore neglect any diagram which does 
not involve a top-quark. This means that the only R-violating couplings which are relevant 
to our discussion are A^3^ (9 parameters) and Ag^^ (3 parameters). 

Furthermore, the values of the top-quark diagrams at m| — > provide an excellent 
approximation to the full integral. Henceforth we work in this approximation. (We note 
that the diagrams carrying only massless fermions would vanish in this limit even if we had 
previously retained them.) 

In the limit m| — >■ 0, the four classes of diagrams can be written in terms of the 1/e pole 
piece and two independent functions of the fermion-scalar mass ratio x = m'j/m'j. which we 
call f{x) and g{x). Their explicit forms are shown in the Appendix. We note that g{x), 
the function which appears in the two-scalar one-fermion, and the fermion wavefunction 
renormalization diagrams, vanishes rapidly as a; —>■ 1. Thus, the finite pieces of the the 
two-scalar one-fermion and the wavefunction renormalization diagrams may be neglected 
for small scalar-fermion mass splittings. Note further that if the poles of these diagrams 
cancel (as is the case for diagrams involving gluinos, for example), the m| = finite pieces 
will also cancel. These considerations apply equally to R-conserving corrections leading to 
significant simplifications in their contributions as well, the details of which will be discussed 
in Sec. 0. 

III. CORRECTIONS FROM THE A' INTERACTIONS 

The R-parity violating A' interactions expressed in terms of the component fields take 
the form 



4 



ijk 



{GihdkRUjL + UjidkRCiL + dlj^el^UjL) 



h.c. 



(3.1) 



As discussed in the previous section, the dominant corrections to hadronic Z decays from 
these interactions are those which involve the top-quark. These are shown in Fig. 0. (The 
diagrams with an internal top-quark and external leptons were considered in Ref. 0.) This 
necessarily means that only the couplings of the right-handed down-type quarks di^ to the 
Z are corrected in our approximation. 

Using notation established in Ref. 0, the corrections to the Z decay amplitude from 
these diagrams are 



I \3k I 



-^ — Z'^ip + q) qkn{p)li^<lkM) 

,2 



COS 9w 

(la) : 2/^6^(724 (0,0,m^;mt, me,^,me^^ 
(16): hu^ ((i- 2)^24 (0,0,m|;me^^,mt,mj 
-m\C23 (O, 0, m|; m-e^^.mt, rut 
(Ic) : -hujim^Co{o,0,ml;mg^^,mt,mt 
(lci) + (le) : 2hd^Bi{0;mt,meJ 



where 



hfL = hf - Qf sin^ e 



w, 



(3.2) 



(3.3) 



The tree level amplitude is /i^.^ times the expression in the square brackets. These corrections 
can be expressed clS cl shift in the coupling hd^^: 



6h. 



i3k 



lA 



/ |2 

i3k I 



+hu^ {{d-2)C2A {me^^,mt,mt) - m\C23 {me^,^,mt,mt)} 



(3.4) 



Henceforth we assume a common slepton mass mg j- 

4 = 

IKskl'H^) 



mg, i = 1,2,3. The full expression 
for 6hi3k is well approximated by the leading m| = piece of the expansion in the Z mass: 

1 



5h 



i3k 



2(47r) 



where 



F{x) = f{x) +g{x) 



X 



1 — X 



1 + 



1 



1 — X 



\nx 



X 



mj' 



(3.5) 



(3.6) 



For mg = lOOGeV, this becomes 
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5h, 



3k 



-0.215% |A^3fc|^ 



(3.7) 



The full shift to the coupling of the quark d^j^ to the Z due to R-violating A' interactions is 
then obtained by summing over the slepton generation index i 



^ -0.215% j:\Ksk\' 



(3.8) 
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Observe that this is a different combination of A' couplings than the combination Y,k l^isk 
which is constrained by lepton universality in Ref. R]. 



IV. CORRECTIONS FROM THE A" INTERACTIONS 

The R-parity violating A" interactions expressed in terms of the component fields take 
the form 



AC' 



2'^ijk 



uld]dl + u'ld*dl + M* c?5(i; 



i ^j^k 



+ h.c. 



(4.1) 



The SU{3) color indices are suppressed. Note that X'-ji^ is antisymmetric in the last two 
indices due to color anti-symmetrization. 

Again, the corrections involving a top-quark are necessarily those with a right-handed 
down-type quark on the external legs as shown in Fig. ^j. Their respective contributions to 
the amplitude are: 



-2 1 A 



// 1 2 

3jk\ 



COS 9 



w 



Z>'{p + q)qjR{p)-f^qjR{q) 



(2a) : -2/1^^(724 (0,0,m|;mt,mj^^,mj^ 



{d - 2)^24 (O, 0, my, mj^^, mt, rrit 
-mlC23 (O, 0, m|; mj^^, rut, rrit 



(2c) : K,m^,Co (o, 



,mz;m^ mt,mt 



(4.2) 



(2d) + (2e) : 2/i,^5i (O; m^, m^^ 
The common leading factor of 2 in this equation is a consequence of the identity 

^abc^a'bc '^^aa'- 

These corrections shift the coupling of the right-handed down-type quark to the Z by 



Sh 



3jk 



-2 1 A' 



// 1 2 

■3jk\ 



■hun {{d- 2)^24 (mj^^, mt, Tilt) - mlC23 (m^^^,mt, m^) } 
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(4.3) 



We assume a common squark mass rnj^^^ = mj, k = 1,2,3. As in the A' case, we have 
neglected all diagrams which vanish in the limit mz — in the above expression. Applying 
the same approximation to the leading diagrams leaves: 

Shs.k ^ j^,\K,,mx) (4.4) 

with X = mf/m'^. For = lOOGeV, this becomes 

6h,k^-0A3%\\'i^,\\ (4.5) 

The full shift to the coupling of the quark rfj^ to the Z due to R-violating A" interactions 
is then obtained by summing over the slepton generation index k 

^ \a?>%Y.\\1^,\\ (4.6) 

k 

In contrast to the A' case. A" interactions do not correct any of the lepton couplings 
to the Z. Thus they do not give rise to lepton universality violations, and no additional 
constraints on A" couplings arise from analysis of the lepton sector. Thus, all significant 
R-violating A" shifts to Z pole observables appear as shifts to the effective coupling of the 
Z to right-handed down-type quarks. 



V. CORRECTIONS FROM R CONSERVING INTERACTIONS 

As stressed in the introduction, in order to isolate the effects of R-violating interactions 
we must properly parametrize the R-conserving radiative corrections (a partial study of 
these effects has also been performed in ||25|). We work in the limit of degenerate sfermion 
masses and tan/3 not large. In this limit, only two parameters are necessary to account for 
R-conserving effects. 

We list all relevant vertex corrections from R-conserving MSSM interactions: 

chargino— sfermion loops : 

The fermion interactions with the gaugino component of the chargino can generate 
substantial corrections to the left-handed couplings of all the fermions (Fig. |^). The 
correction to the up-type and down-type quark couplings from the diagrams shown 
in Figs. PD,c,d are proportional to 

^^g.c.d) ^ 2/i„^ C'24("^x'"^"l>"^«l) + hd^Bi{m^,miiJ, (5.1) 

where the dependence on the external momenta have been suppressed. In the limit 
m| — i> and = , it is clear (see Appendix) that 
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Sh(^b,c4) ^ .^/.g'A^) 0, (/.^^ _h,J = il- sin' Ow). (5.2) 

A similar relation exists for the correction to the leptonic vertices provided m^^ = mg^ . 
In addition, the diagram of Fig. ^ changes sign with the isospin of the final-state 
fermion. As a result, the combined contribution from all the diagrams in Fig. |^ is 
proportional to the isospin of the final-state fermion but otherwise universal in the 
limit that all the (left-handed) squark and slepton masses are degenerate. 

A shift proportional to the isospin can be written as an overall multiplicative change 
in the coupling and a shift in the effective value of sin^ 9w'- 

hf^ = l3f{l + 6) -Qsin'Ow = 

hf^= -Q sin' Ow = {1 + 6)1^ -Q "-^^ff) ^^-^^ 

Since we utilize only observables which are ratios of couplings, the multiplicative cor- 
rection cancels and only the shift in sin^ 6w is measurable. 

charged Higgs— top (Higgsino— stop) loops : 

The charged Higgs-top corrections (Figs. ^,d,e) and the supersymmetrized versions 
of these diagrams containing the chargino-right handed stop loops (Figs. Hb,c), as well 
as the corresponding wavefunction renormalizations are peculiar to the bi final states. 
Higgs and higgsino couplings to the and other fermion final states are suppressed by 
the light fermion masses and are neglected. The model generically predicts Sh^^^^'^ > 0. 
This is a result of the fact that (1) the charged Higgs contribution is always positive 
and (2) explicit calculation shows that the leading p' /m? Higgsino contribution 




vanishes.]^ 

gluino— squark loops : 

The leading (m| = 0) contributions from the diagrams shown in Fig. ^ vanish as a 
result of the relation (see Appendix) 



2 (724(0, 0, m|; mg, mg, niq) + Bi (0; m^, rriq 



0. (5.4) 



Even for the subleading terms, in the limit of degenerate squark masses the gluino- 
squark loops induce only a universal shift to all of the Z-quark couplings. Just as for 
QCD corrections, this shift cancels in the ratios of hadronic widths and asymmetries 
and therefore does not enter into our analysis. 

neutralino— sfermion loops : 

The gaugino component of the neutralino generates shifts to all of the couplings 



^For the purpose of this analysis, we treat Higgsino and gaugino contributions separately. 
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(Fig. 1^). However, the only (potentially) significant diagrams are all of the two scalar, 
three sparticle variety. When these are combined with wavefunction renormalization 
diagrams, the leading m| = pieces cancel in the sum (see Eq. |5.4| ). Therefore, these 
can be neglected altogether. 

oblique corrections : 

In addition to all these vertex corrections, R-conserving interactions can also affect Z- 
peak observables through vacuum polarization diagrams, aka the oblique corrections. 
Oblique corrections can all be subsumed into a shift in the p parameter and the effective 
value of sin^6'vi/, the first of which cancels in all of the observables that we consider 

Thus, in our approximation, the only R-conserving effects we need to consider are (1) a shift 
in the left-handed h coupling from charged Higgs/Higgsino corrections, and (2) a universal 
shift in the effective value of sin^ dy^ which subsumes the isospin-proportional correction due 
to chargino-squark loops as well as the oblique corrections. 



VI. FIT TO THE DATA 

As we have seen, the R-violating A' couplings correct the left-handed couplings of the 
charged leptons and the right-handed couplings of the down-type quarks while the A" cou- 
plings correct the right-handed couplings of the down-type quarks only. The R-conserving 
chargino-sfermion correction is absorbed into a universal shift of sin^ Qy^ (provided that all 
the sfermions are degenerate) while the Higgs-top correction is only relevant for the left- 
handed h quark. Since we have already discussed the limits placed on the A' couplings from 
the leptonic observables in a previous paper 0, we will concentrate on the corrections to 
the quark observables and perform a fit which encompassed both the A' and A" cases. 

In order to constrain the size of these corrections we will use the ratios of the hadronic 
parital widths 

c,6) 

(g = u, d, s) 



and the parity-violating asymmetry parameters 

-h^ 3 

These observables have the convenient property that (1) they are insensitive to QCD and 
gluino-squark corrections, and (2) the only dependence on oblique corrections (vacuum 
polarizations) comes from a shift in the effective value of sin^ 6w This will permit us to 



Rr, 



R' 



had 



K + Kn 



--u,d,s,c,b 



K +K 



E 

q=u,dy 



{h 
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constrain the parameters we are interested in without comphcating the fit procedure by 
introducing gluon/gluino corrections or corrections to the p parameter. 

Of the leptonic observables, we will include the ratio of electron to neutrino widths 
Ry/e = T^p/T^+e- and the electron asymmetry parameters and Ap-B{e) = ^A^ to help 
constrain the universal R-conserving and oblique corrections. Corrections to Ae must be 
considered in any case since it is present in the hadronic observables AfbIq')- Though the 
left-handed lepton couplings receive corrections from the A' interactions, the size of the 
correction particular to the electron is already so tightly constrained to be small by other 
experiments that we can neglect it entirely. We drop all /i or r dependent observables from 
our fit so that we can use the result to constrain both the A' and A" cases. 

In table H we list the experimental data we use in our fit with correlation matrices shown 
in tables and We caution the reader that many of these numbers are preliminary 
results announced during the summer 1999 conferences so they, and our resulting fit derived 
from them, may be subject to change. Some comments are in order: 

1. The ratio 

r - h'^ 

J- e+e "-eL ^ '^e^ 



is calculated from the first six Z-lineshape observables. Its correlation to AFB(e) is 
+28%. Its correlations to the \i and r observables, which we drop, are negligibly small. 

2. The r polarization data has been updated from Ref. |l^ with new numbers from 
DELPHI [0,13. We keep only A^ and drop Ar. 

3. The SLD value of Alr (which is the same thing as A^ is from hadronic events only. A^. 
is from the leptonic events. Its correlations to the dropped A^ and Ar are negligibly 
weak. (The errors are dominated by statistics [|16|-) 

4. The OPAL measurements of i?^, Afb(s), and Ap^{u) assume Standard Model values 
of Bl^ = 0.359 and ApBid) = 0.100. To account for the shifts in the down observ- 
ables in our model, the data should be interpreted as constraining the following linear 
combinations: 



R'* = R'^ + 1.83 [ R'a - 0.359 ] 
A*p^{s) = Afb(s) - 0.32 [ AFBid) - 0.100 ] 
A;b(m) = Afb(m) - 1.42 [ AMd) - 0.100 ] 



There is a +31% correlation between Ap-^ls) and A^ 



^FB 



.u 



5. The DELPHI measurement of Apb{s) assumes standard model values of Apsiu) = 
0.0736 and ApBid) = 0.1031. It should be interpreted as a measurement of the follow- 



ing linear combination |T8 

A**b(s) = Afb(s 



0.156 [ Afb{u) - 0.0736 ] - 0.117 [ ApB^d) - 0.1031 
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6. The SLD measurement of Ag [|T9| assumes standard model values for Au, A^ 



To 



and R'^. (The dependence on the heavy flavor variables are weak and negligible, 
account for shifts in these input parameters the measurement should be intepreted as 



constraining pO| 



Al = As- 0.0602 [ Au - 0.668 ] - 0.0467 [ Ad - 0.936 
- 1.32 [ R'^ - 0.280 ] - 1.20 [ R'^ - 0.360 ] 



The heavy flavor data is the combined fit to the LEP and SLD data compiled by the 

±c are shifted 



LEP Electro weak Working Group The central values of A}, and A^ 

compared to the original SLD values of 



Ab = 0.905 ± 0.026 
Ar = 0.634 ±0.027 



Using these numbers instead of those shown in Table | will result in a slightly tighter 
constraint on the R-violating couplings, but we will present the results using the 
LEPEWWG numbers to be on the conservative side. 

We denote the shift in sin^ 9w due to oblique and chargino-sfermion corrections by 6s', 
and the shift from the Higgs interactions specific to the left-handed coupling of the b by 
6h^^^^^. The R-violating shifts specific to the right-handed couplings of the d, s, and b 
quarks are denoted 6hf^, 6hf^, and Shf^. Then the shifts in the couplings of the quarks, the 
electron, and the neutrino are given by: 



= 
= 




5s' 
6s' 




Shu^ = 






= 






Shd^ = 






5hdR = 


hs' 
3 




6hc^ = 


3 




Shcu = 


3 




6hs^ = 






= 


3 




6hb^ = 


3 




Shb^ = 
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The dependence of the observables on these fit parameters can be calculated in a straight- 
forward manner. For instance, we find: 



Ru/e 



h 



hi, + hi 



2he^ + 2h. 



\ "-EL ^ "'en 

0.64 5s^ 



6s' 



or 



5R^/e = 1.1755^ 

where the coefficient has been calculated assuming sin^ 9w 

6A. 



0.2315. Similarly, 



6R': 



-7.61 Ss^ 
-1.63(5.2 

0.151 + 0.242 5/1^ 



- 0.0058 

d-R sn 



^FB 



U 



6A;^{s) 

6A;Us) 
6A: 

SRb 
SRc 

SAMb) 
SAfb{c) 
SA, 
6A. 



3.74 5.2 + 0.262 5/1 
-3.72 5.2 + 0.0558 - 0.174 
-4.15 5.2 + 0.020 - 0.174 
-0.321 5.2 - 0.0444 5/if^ 

0.0392 5.2 -w.wo^yu^./^^^ 
-0.0605 5.2 



1.37 5hf^ 



bR 



0.0396 6h^^ - 0.0396 6hf^ + 0.141 Sh^ 
0.0316 5/if - 0.0316 5/if 

<1r Sfi 

-5.40 5.2 - 0.172 5/if^ - 0.0315 5/i"'^*^' 
-4.175.2 
-0.636 5.2 
-3.45 5.2 



0.771 5/i^f^^ 



0.0316 5/if^ + 0.173 5/i™^s^ 



1.615/1 



bR 



0.295 5/iff^^ 



(6.1) 



Fitting these expressions to the table | data, we obtain: 



5.2 
Shi 
Sht 



6h 



5/iff^ 



bR 



-0.00092 ± 0.00022 
0.081 ±0.077 
0.055 ± 0.043 
0.026 ±0.010 

-0.0031 ±0.0042 



(6.2) 



with the correlation matrix shown in table |^ The quality of the fit was x = 12.0/(16 — 5). 



usmg rrit 



The standard model predictions were obtained using ZFITTER v. 6. 21 [^| 
174.3 GeV |22| and rrih = 300 GeV. Except for 5.2, the best-fit values and uncertainties of the 



parameters are virtually unchanged when the Standard Model Higgs mass is varied between 
100 GeV and 1 TeV. By far the largest contribution to the is from those observables 
(i?i//e, Afb(c) and Ac contribute a combined 8.6) which serve only to compete with Alr in 
determining 5.2. 

In Figs. 1^ through |l^, we show the limits placed on the five parameters by various 
observables projected onto two dimensional planes. Figs. |^, [TT], [12], and |T3| show that the 
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most stringent constraint on Sh^^ comes from A*p^{u), while Figs. H, 0, |14|, and |T5| show 



that the 6hf^ is constrained by A* and Figs. |Tg, |T3|, 0, and |T6| show that dh^l^^"" is 



largely fixed by Rb- It is clear from figs. [T^, |I4|, and |16]that the strongest constraint on Sh^'^ 
comes from Ai, and ApBib). However, a careful look at Fig. |^ shows that the limit on Shf^ 
is strongly correlated with the value of (5s^. Because and other measurements prefer a 
slightly negative 5s^, the preferred value of Shf^ from ApBib) is shifted to the positive side 



VII. LIMITS ON A' AND A" 



Using Eq. [3^71 , we can translate our fit results in Eq. to limits on the A' couplings 
constants: 

^|A:3i|2 = -38 ±36 
= -26 ±20 

^|A:33p = -12.1± 4.7. (7.1) 

i 

The correlations between the fit values of the couplings are relatively small (see table [TVp . 
The Icr [2a) [3a] upper bounds are then 

E|A:3,P<-2 (34) [70] 
E|A:32p<-6 (14) [34] 

E|A:33I'<-7.4(-2.8)[1.9] (7.2) 

i 

This imposes the following {2a) [3a] limits on the individual couplings in the sum: 

\X[,,\ < (5.8) [8.4] 
|A:32| < (3.8) [5.9] 

|A:,33I<( )[1-4] (7.3) 

For i = 1 and i = 2, stronger constraints at the 2a level on the relevant couplings are available 
from other types of experiments ||^, so these constraints fail to improve previous results. The 
strongest constraint is on i = 3, where the best-fit value of the sum of squared couplings 
is negative even at 2a. This constitutes a significant improvement in the upper bound on 
IKssl o^^^ previous bounds on these couplings which were nonzero at the 2a level. These 
results are complementary to those obtained in Ref. [Q, in which a different combination of 
A' couplings was constrained. In particular, for the A'33 couplings, the a from the lepton 
universality constraints is much smaller, but the best-fit value of the squared couplings from 
the hadronic constraint is negative by an even greater statistical significance. Q 



strong independent constraint on A-gg will be available from the measurement of the invisible 
width of T resonance 
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Next, we consider the constraints on the A" couphngs. These couphngs have hitherto 
been constrained by experiment only weakly or not at all. Using Eq. |4.5|, we can translate 
Eq. into the bounds: 



5:|A'3;,p = -19±18 
El^m.r = -13±10 



X:|A'3'3,r = -6.0±2.3. (7.4) 

k 

Again, the correlations between these constraints are relatively weak, so we neglect them 
henceforth. The la {2a) [3a] upper bounds are then: 

ElAmr<-l ( 17) [35] 

k 



El^mP<-3 ( 7) [17] 

k 

El^mP<-3.7(-1.4)[0.9]. (7.5) 

k 

Recall that the A" couplings are antisymmetric in the last two indices; thus, each of the 
sums above consists of only two terms. The {2a) [3a] upper bounds on the individual A" 
couplings are then: 

lA^'^il < (2.7) [4.1] 

|A'33J<( )[0.96]. (7.6) 

Thus, A33^ and A332 are excluded at the 2a level, and A321 is excluded at the la level. These 
bounds significantly improve the la bound of |A33^| < 0.50 from Ri [§|,^. 

These improvements on the bounds of A' and A" are a consequence of the fact that while 
the data prefers a positive shift in the right-handed coupling of the b, which is non-zero by 
2.6a, both A' and A" corrections shift the coupling in the negative direction. This situation 



is mitigated neither by introducing sfermion mass splittings nor by increasing tan/? [26 . 



VIII. BAYESIAN CONFIDENCE INTERVALS FOR A' AND A" 

In the previous section we performed a classical statistical analysis, i.e. we performed 
a fit to the data without any a priori assumptions about the viability of the model. In 
particular, we made no assumptions about the signs of the coupling shifts when fitting the 
data. As a consequence, the best-fit values for the squares of the R-violating couplings were 
negative, resulting in strong la and 2a bounds. 

An alternate method for calculating confidence levels is to use Bayesian statistical anal- 
ysis. This technique assumes that R-violating SUSY is the correct underlying theory, and 
therefore that the shifts to the right-handed couplings are only permitted to be negative and 
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Sh^^^^^ positive.[] The resulting confidence intervals for the couplings squared are positive, 
and the preferred values are those of the Standard Model (i.e. zero). 

However, care should be taken when using these bounds, since they hide the fact that 
the of thfi corresponding fit is quite large even at low confidence levels. The probability 
of these bounds arising as a result of statistical fluctuations is therefore quite small. 

Below we list the 68% (95%) confldence levels from the constrained fit: 



ibn ^ -0.0046 (-0.010) 
5, > -0.031 (-0.064) 

> -0.061 (-0.123); (8.1) 

The corresponding confidence limits on the couplings are: 

IA33J < 1-0 (1.5) 

IA'3'211 <2.7(3.9) (8.2) 

|A^33| < 1.4 (2.2) 
|A-32l < 3.8 (5.6) 

|A:3i|<5.2(7.6); (8.3) 

The best-fit value for 6h^^^^^ is negative. However, the model generically predicts a 
positive The best-fit value of 6h^^^^'' consistent with the model is zero; as a result 

of this tension, the corresponding increases even further. 

To be quantitative concerning the large of the constrained fit confidence intervals, 
we present the following example. The corresponding to the 68% and 95% confidence 
intervals for Shf^ are: 



68% : X /DOF = 26.1/(16 - 1) ^ probability = 3.7% 







95% : X /DOF = 33.6/(16 - 1) ^ probability = 0.4% (8.4) 

We see explicitly that the bounds obtained using the Bayesian analysis are weak, but 
the associated with these bounds is uncomfortably large. If the error bars continue to 
shrink and the central values stay unchanged, the relevance of the constrained fit bounds 
must be questioned. 



IX. SUMMARY AND CONCLUSIONS 

We find that the hadronic Z-decay data from LEP and SLD can be used to place 
significant constraints on the size of R-parity violating A' and A" couplings. This is possible 
because the dominant R-violating interactions correct the couplings of the right-handed 



^We take the prior probability for the coupling shifts to be uniform on the region permitted by 
the theory and zero elsewhere. 
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down-type quarks only while the dominant R-conserving MSSM interactions correct only 
the left-handed couplings. The parity violating asymmetry parameters Ag are particularly 
sensitive to shifts in the right-handed quark couplings while blind to shifts in the left-handed 
couplings. This allows us to constrain the R-violating interactions independently from the 
R-conserving sector. 

Current data prefer a shift in right-handed quark couplings opposite to the direction 
predicted by the theory. As a consequence, all of the R-violating shifts considered in this 
work are excluded at the la level. In the A' case the strongest bound is on the A -33, which 
are excluded at 2a and on which we have set the 3a bound 

|A:33| < 1.4. (9.1) 

For the A" case, the A33^ and A332 couplings are excluded at the 2a level, and A32j^ is excluded 
at la. The {2a) [3a] upper bounds are 

IA:;^,! < (2.7) [4.1] 

|A'3'3J < ( ) [0.96]. (9.2) 

All bounds are calculated assuming a common sfermion mass of 100 GeV. For larger (com- 
mon) sfermion masses the above bounds may be interpreted as bounds on (|A'|, |A"|) x 
F{x)/ F{xo), where F{x) is defined in Eq. and Xq = (Yoc^tvp'- have also performed 
a Bayesian statistical analysis and obtained corresponding confidence levels. 

Generically, R-violating interactions reduce the magnitude of the couplings of the the 
right-handed quarks to the Z and leave the left-handed couplings unchanged. Current 
LEP/SLD data prefers shifts which increase the magnitude of the right handed coupling, 
to the extent that even the Standard Model prediction is in only marginal agreement with 
the data. Future reductions in the experimental uncertainties in the asymmetry parameters 
without changes in the central values would eventually rule out both the standard model 
and the MSSM with R-violating couplings of either the A' or A" variety. 
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APPENDIX: FEYNMAN INTEGRALS 



The integrals we use here are defined exphcitly in In the approximation m| = 0, the 
one-loop diagrams which appear in this work are proportional to the following expressions: 




{d - 2) C2A (0, 0, m|; ms,mf,mf^ - m| Caa (o, 0, m|; ms^mf^m/^ 



{Any 



1 / TTl 



+ f{x) 



(Al) 



where 



oc 2 C24 (0, 0, 171^', rrif, rris, m, 



(4vr) 



1 / TTl 

2 ^'-'"7/, 



(A2) 




oc mj Co (0, 0, m|; ms,mf, rrif 



(An)' 



[f{x)+g{x) 



(A3) 




oc -Bi (0; nif, 



(47r)^ 



1 / TTl 



(A4) 



fix) 
9{x) 



4(1 

-mx -\ — ; — - 

2 4(1 



X — 1 — 2 In a; 



-(l-a;)(l-3x) + 2x^ Inx 



for X = mi /ml. For x — > 1 (degenerate scalar and fermion masses), 



fix) 
9{x) 



1 x-1 



2 6 

X — 1 



+ 



(A5) 



(A6) 



For X — > (the decoupling limit of heavy scalar masses) , 

/(x) ^ ^Inx + ^H 

^(x) ^-ilnx-i + --- ~ -f{x) 



(A7) 



17 



TABLES 



Observable 



Reference 



Measured Value 



ZFITTER Prediction 



Z lineshape variables 






|14| 

121 

|14| 

121 

|14| 





91.1872 ±0.0021 GeV 
2.4944 ± 0.0024 GeV 
41.544 ± 0.037 nb 
20.803 ± 0.049 
20.786 ± 0.033 
20.764 ± 0.045 
0.0145 ± 0.0024 
0.0167 ±0.0013 
0.0188 ± 0.0017 
1.9755 ±0.0080 



input 



0.0152 



1.9916 



T polarization at LEP 











0.1483 ±0.0051 
0.1424 ±0.0044 



0.1423 



SLD left-right asymmetries 




0.15108 ±0.00218 
0.1558 ±0.0064 
0.137 ±0.016 
0.142 ±0.016 



0.1423 
0.1423 



light quark flavor 



A 



R'; [OPAL] 
[OPAL] 



FBI 



A*p^{u) [OPAL] 
[DELPHI] 
At [SLD] 



^FB 





[17| 

121 

m 

m 



0.392 ± 0.062 
0.075 ± 0.029 
0.086 ± 0.037 
0.1008 ±0.0120 
0.85 ± 0.092 



0.360 
0.100 
0.071 
0.1006 
0.935 



heavy quark flavor 






121 






in 



0.21642 ± 0.00073 
0.1674 ±0.0038 
0.0988 ± 0.0020 
0.0692 ± 0.0037 
0.911 ±0.025 
0.630 ± 0.026 



0.21583 
0.1722 
0.0997 
0.0711 
0.934 
0.666 



TABLE I. LEP/SLD observables and their Standard Model predictions. The ratio 

The Standard Model pre- 
GeV H, ma = 300 GeV, 



R,y/e = ^up/^e+e- '^^^ Calculated from the Z-lineshape observables. 
dictions were calculated using ZFITTER v.6.21 |2ll with mt = 174.3 
and as{mz) = 0.120 as input. 
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mz Tz 




"^had 


He 






^FB(e) 




A / ^ 

^fb(t) 


mz 


1.000 -0.008 


-0.050 


0.073 


0.001 


0.002 


-0.015 


0.046 


0.034 




1.000 


—0.284 


-0.006 


0.008 


0.000 


-0.002 


0.002 


—0.003 




'''had 




1.000 


0.109 


0.137 


0.100 


0.008 


0.001 


0.007 


Re 






1.000 


0.070 


0.044 


-0.356 


0.023 


0.016 


Rfj, 








1.000 


0.072 


0.005 


0.006 


0.004 


Rr 










1.000 


0.003 


-0.003 


0.010 


^FB(e) 












1.000 


-0.026 


-0.020 


^fb(/^) 














1.000 


0.045 


^fb(t) 
















1.000 



TABLE II. The correlation of the Z hneshape variables at LEP. 





Rb 


Rc 


AMb) 


Afb(c) 


Ab 


Ac 


Rb 


1.00 


-0.14 


-0.03 


0.01 


-0.03 


0.02 


Rc 




1.00 


0.05 


-0.05 


0.02 


-0.02 


AMb) 






1.00 


0.09 


0.02 


0.00 


Afb{c) 








1.00 


-0.01 


0.03 


Ab 










1.00 


0.15 


Ac 












1.00 



TABLE III. The correlation of the heavy flavor variables from LEP/SLD. 







6hf 




5h^ 

Or 


br, 




1.00 


0.01 


-0.06 


-0.42 


-0.15 


< 




1.00 


-0.30 


0.09 


-0.75 






1.00 


0.05 


-0.22 










1.00 


0.30 
1.00 



TABLE IV. The correlation matrix of the fit parameters. 
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FIGURES 





(d) (e) 
FIG. 1. One-loop corrections to Z — > dfe^d/j^ involving the R-parity violating A' couplings. 
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(c) (d) 

FIG. 3. Examples of leading IPI R-parity conserving chargino-sfermion contributions sub- 
sumed into (5s^. There are analogous diagrams with ul quark final states and with lepton final 
states. 
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(a) 



(b) 



FIG. 5. Gluino-squaxk corrections to the Zqq vertex. 




(a) (b) 

FIG. 6. Neutralino-squark corrections to the Zqq vertex. 
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0.5 



Ae(SLD) 



Sh^^ 0.0 K 



-0.5 — . 



-1.0 




-0,01 -0.005 



FIG. 7. Constraints in the Ss'^-Sh^^ plane from various observables 



1.0 



0.5 



SR 



0.0 



-0.5 



-1.0 



-0,01 




0.01 



FIG. 8. Constraints in the 5s^-5hf^ plane from various observables 
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0.10 



0.05 



5hf^ 0.00 ^= 



-0.05 — 



-0.10 




-0.01 -0.005 



FIG. 9. Constraints in the Ss^-dhf^ plane from various observables 



0.10 



0.05 — 



0.00 



-0.05 



-0.10 




FIG. 10. Constraints in the Ss^-Sh^^^^^ plane from various observables. 
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1.0 



0.5 



Shf^ 0.0 



-0.5 — 



-1.0 




FIG. 11. Constraints in the 5h^^-5hf^ plane from various observables. 
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0.05 — 



0.00 
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-0.10 




FIG. 12. Constraints in the 5h^ ^^h? plane from various observables. 
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0.10 



0.05 



shf:^' 0.00 - 



-0.05 — 



-0.10 




FIG. 13. Constraints in the 5h^^-5hf^ plane from various observables. 
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0.05 — 



0.00 



-0.05 



-0.10 




5hl 



FIG. 14. Constraints in the 5h^ ^^h? plane from various observables. 
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FIG. 16. Constraints in the Shf^-Sh^^^^ plane from various observables. 
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